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Abstract Numerical simulations of the helical {m = 1) kink instability of an 
^ \ arched, line-tied flux rope demonstrate that the helical deformation enforces 

reconnection between the legs of the rope if modes with two helical turns are 
CO \ dominant as a result of high initial twist in the range ^ > Gtt. Such reconnection 

is complex, involving also the ambient field. In addition to breaking up the 
original rope, it can form a new, low-lying, less twisted fiux rope. The new fiux 
q/ ■ rope is pushed downward by the reconnection outfiow, which typically forces it 

^Q \ to break as well by reconnecting with the ambient field. The top part of the 

original rope, largely rooted in the sources of the ambient fiux after the break- 
up, can fully erupt or be halted at low heights, producing a "failed eruption." 
The helical current sheet associated with the instability is squeezed between 
the approaching legs, temporarily forming a double current sheet. The leg- leg 
reconnection proceeds at a high rate, producing sufficiently strong electric fields 
^ I that it would be able to accelerate particles. It may also form plasmoids, or 

plasmoid-like structures, which trap energetic particles and propagate out of the 
reconnection region up to the top of the erupting fiux rope along the helical 
^ . current sheet. The kinking of a highly twisted fiux rope involving leg- leg recon- 

f^ I nection can explain key features of an eruptive but partially occulted solar fiare 

"^ . on 18 April 2001, which ejected a relatively compact hard X-ray and microwave 

^^ I source and was associated with a fast coronal mass ejection. 
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1. Introduction 



When a line-tied, arched magnetic flux rope becomes unstable with respect to 
the helical kink mode, i.e., to displacements with azimuthal wavenumber tti = 1, 
several effects occur which are relevant to the understanding of solar eruptions 
dSakurai, 19761 |Fan ^^^ Gibson, 20031 [Torok, Kliem, and Titov, 2004D . First, the| 
apex of the rope will ascend if the initial perturbation points upward, which is 
usually the case for solar eruptions due to a slow rise of the structure prior to the 
main acceleration. Second, the ensuing fast rise is initially exponential or a power 
law close to an exponential, as observed to be typical for fast coronal mass ejec- 
tions (CMEs) ( jVrsnak, 20QT1 [Schrijver et a/., 2008a| . Third, the rope wih writhe 
by the conversion of twist helicity into writhe helicity. This leads to a rotation of 
the apex about the direction of ascent and forms a characteristic inverse gamma 
shape or O shape (depending on the perspective of the observer). Such shapes 
are often observed and have been taken as evidence for the occurrence of the 
instability (e.a. JRustll2003|IJi et a/.ll2003| Romano, Contarino, and Zuccarello 



[20031 iRust and LaBontel [20051 [Torok and Klieml [20051 iFanl [20051 [Zhou et al. 



[20061 [Gilbert, Alexander, andTiul [20071 ICho et aU [2009). Fourth, two current 
sheets will be formed and exponentially steepened, so that the onset of magnetic 
reconnection is likely: a vertical (flare) current sheet under the rope resulting 
from the rope's rise; and a helical current sheet, wrapped around the rope and 
passing over the rope's apex, created as a consequence of the helical displace- 
ment. Reconnection in these current sheets influences the eruption strongly. 
It leads to a variety of different dynamical behaviors, including eject ive and 
conflned eruptions ( Torok and Kliem, 2005^ as well as the break-up and refor- 
mation of the rope, which we will address in the present paper. Fifth, sigmoids 
are likely to form, since the fleld lines that thread either of the current sheets 
are sigmoidal in projection (Titov and Demoulin, 1999 Fan and Gibson, 2003 



Kliem, Titov, and Torok, 2004| [Gibson et al, 2004[). The orientation of the sig- 



moidal fleld lines underneath the flux rope, where sigmoids are known to form, 
correlates with the chirality of the fleld in agreement with the observations 



( Rust and Kumar, 1996 ; Pevtsov, Canfleld and McClymont, 1997 Green et a/., 2007| ).| 

The literature reports a number of events that suggest the occurrence of the 
helical kink instability through observational evidence for both the existence of 
considerable twist (exceeding one fleld line turn about the axis of the structure) 
and the development of writhe in the course of the eruption. Examples are the 



events of 19 July 2000 ( Romano, Contarin o, and Zuccarello, 2003), 18 August 
1980 dVrsnak et a/., 1993[), and 12 June 200 3 (Liu and Alexander^OOSj) ; see also 
the systematic study by 'Vrsn ak, Ruzdjak,^and Rompolt| (;1991j) . 

In spite of such considerable principal agreement and support, it is possi- 
ble that the instability actually occurs less frequently than indicated by the 
above observations, because there exists a second physical process that can 
cause the writhing of a rising flux rope into helical shape. As has been pointed 
out by [Isenberg and Forbes] (|2007p , the shear fleld component of the ambient 
fleld (pointing along the photospheric polarity inversion line), crossed with the 
vertical current component in the legs of a flux rope which is displaced upwards 
out of equilibrium, gives a sideways pointing Lorentz force on the flux rope 
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legs. This force yields a writhing in the same direction as the kink- mode-driven 
writhing for any process that Hfts the rope out of equilibrium. Basically, the 
same observational consequences as listed above can be expected, including the 
reduction of any possibly existing initial twist into the developing writhe, due 
to helicity conservation. 

Therefore, to evaluate the relevance of the helical kink instability for solar 
eruptions, it is important to distinguish between the two possible causes of the 
writhing into helical shape. The direct approach — computing the coronal field by 
extrapolation of vector magnetograms and comparing the obtained twist to the 
threshold of the helical kink mode — is still hardly practicable because present 
vector magnetograms and nonlinear force- free extrapolation codes do not appear 
to have reached the degree of consistency required for a reliable computation 
(|Metcalf et al.. 2QQ8,,Schrijver et al.^ 2QQ8b) , and because a comprehensive para- 
metric study of the instability threshold has not yet been completed (although 
some knowledge is already available, see'T orok, Kliem, and Titov[ 120041 and the 
references therein). Hence, it is of interest to identify features that permit an 
observational discrimination. 

The writhing by the shear field component has the property that the sideways 
pointing Lorentz force acting on the flux rope legs has the same sign all along one 
leg, and the opposite sign along the other leg; it changes sign only at the apex. 
This implies that the displaced flux rope cannot develop more than one helical 
turn. In terms of a normalized axial wavenumber, k' = kl/2'K = //A, where / is 
the length of the flux rope and A is the wavelength of the perturbation, this can 
be expressed as /c' < 1. 

In contrast, the helical kink instability can have any wavenumber in principle. 
Its dominant wavenumber is set by the twist ^ = IB^jy/rBz = 27rA/', where N is 
the number of fleld line turns in the flux rope, in general averaged over its cross 
section. Stability analyses indicate that maximum growth occurs in the range 
A ~ (l-2)/p, where /p = l/N is the pitch length, the axial length required for 
one fleld line turn ( Linton, Longcope, and Fisher, 19"96{ [Linton et al.^ 1998[ ). For 



a force-free flux rope with uniform radial twist profile, for example, maximum 
growth occurs for A = 2/p. In this case, the k' r^ 1 and k' ^ 2 modes are dominant 
for $ ~ 47r and ^ ~ Stt, respectively, so that the transition between them occurs 
in this range of twist values. Such high values render the occurrence of /c' > 1 
modes significantly less likely than the occurrence of single helical modes. For 
twist profiles that lead to maximum growth closer to A = /p, the transition will 
occur between twists of 27r and 47r, i.e., closer to the kink instability threshold, 
so it could be more common. However, it is not clear whether such nonuniform 
profiles are often formed in the solar atmosphere. The observations of erupting 
filaments rather tend to indicate A:' < 1 helical modes in the majority of events. 
If modes k^ > 1 are observed, the helical kink instability is implied. Projected 
onto the plane of the sky, this manifests itself as the occurrence of more than 
one infiection point of the kinked fiux tube. For events near Sun center, an S 
shape with an overturn at one or both ends results (see the kinking fiux rope 
in Figure [2] below at t = 11 for an illustration). In fortunate cases the effect 
is apparent even in limb events, for example in the failed filament eruption 
of 27 May 2002, which showed an infiection point in the middle of the front 
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leg (see the middle panels of Fig. 1 in ITorok and Kliem) l2QQ5p . If k' is clearly 
larger than unity, then the multiple helical shape is relatively clearly indicated 
for any perspective; an example is the 19 July 2000 filament eruption analyzed 



in ^Romano, Contarino, and Zuccarello] (j2003p . Such direct observations of the 



helical mode with k' > 1 are not very numerous, however. 

Another possible consequence oi k' > 1 mode dominance is an approach 
of the rising flux rope's legs so that they can come into contact and begin to 
reconnect. Such an evolution is not triggered by a /c' < 1 helical displacement, 
which rather causes the legs of the flux rope to move apart. It is expected to 
have distinct observational consequences. The reconnect ion can cause the break- 
up of the original flux rope as suggested bv lCho et a/.l(|2009p . Furthermore, the 
simulations in Section [2] demonstrate that such break-up can lead to the refor- 
mation of a (less twisted) flux rope from the two bottom parts of the original 
rope, which facilitates the quick reformation of a filament after an eruption. 
Finally, the energetic particles produced by the reconnection can have unique 
radiation signatures if the helical current sheet is involved. This was found in 
the companion paper to the present study ( [Karlicky and Kliem] 12010) hereafter 
Paper I), which analyzes the observations of an eruptive flare on 18 April 2001 
in NOAA Active Region (AR) 9415. This partly occulted eruption formed a 
rapidly ascending microwave source of strongly helical (inverse gamma) shape 
(shown also in Figure [7] below) and presents evidence for the acceleration of 
nonthermal particles at or near the crossing point of the flux loop legs (as seen 
in the projection onto the plane of the sky). Moreover, the particles were effi- 
ciently trapped in superimposed moving compact sources, which are suggestive 
of plasmoids formed by reconnection between the legs and moving in the helical 
current sheet. 

In the present paper, we demonstrate the reconnection between the legs of a 
kinking flux rope using magnetohydrodynamic (MHD) simulations. We estimate 
the minimum twist enabling this process and study its dynamics. The simulation 
results are then discussed in relation to the observations of the eruptive flare on 
18 April 2001 which were found in Paper I to be suggestive of such reconnection. 



2. Helical Kink Instability of a Highly Twisted Flux Rope 

The zero-beta, ideal MHD simulations presented in this paper are all largely 
similar to those in ITorok and Klieml(j2005p . They use the force- free equilibrium 
of a toroidal flux rope by ITitov and DemoulinI (|1999p as initial condition and 
employ similar numerical parameters (a Cartesian box of size 10 x 10 x 20, 
discretized by a stretched grid of resolution A = 0.03 in the central part, and 
very similar small viscosity and weak spatial smoothing, as required to maintain 
numerical stability). Magnetic reconnection occurs in these simulations as a 
consequence of numerical diffusion where current layer (s) steepen in response 
to the development of ideal MHD instability. 

Most geometrical parameters of the rope, as well as the normalization, are 
also set as in these previous simulations: major torus radius R = 1.83, depth of 
torus center d = 0.83, distance of the sources of the external poloidal field (the 
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Figure 1. Flux rope shape for a range of initial twist values and uniform apex height of twice 
the initial apex height (reached at t = 34, 19, 13, and 12rA in the simulations). The volume 
1 X 4 X 4 is shown; the stretching in the x direction helps to visualize the helical shapes. All 
fleld lines are started from circles of radius 0.2a centered at the footpoints of the flux rope. 
The orthogonal projections show the magnetic axis of the rope. 



field components perpendicular to the fiux rope axis) L = 0.83. The line current 
in the model, which produces the external toroidal field, is reduced by a factor 
20 in comparison to the simulation of a failed eruption in that previous paper, 
to enable eruptive behavior ( [Roussev et a/., 2003[ [Torok and Kliem, 2005D . The 
toroidal field component points along the fiux rope axis and is essentially also 
the shear field component, since the polarity inversion line of the field in the 
bottom plane runs approximately parallel to the flux rope under its apex. The 
external field is thus mainly poloidal, with the ratio between toroidal and poloidal 
components at the fiux rope apex being Bq^/Bq^ = 0.075. Consequently, any 
development of writhe is nearly exclusively due to the helical kink mode. Gravity 
is neglected, and the initial density is derived from the initial magnetic field 
by Po = |Bo(x)p/^, to model a corona with an Alfven velocity that decreases 
slowly with distance from the flux concentrations. Initially, the flux rope lies in 
the plane {x = 0}, and its field strength is normalized to unity at the apex, 
^0 = |B(0, 0, 1, 0)1 = 1. This gives an initial Alfven velocity Vao = ^oPo — ^ 
at the apex, which serves as the velocity unit and defines the time unit through 

TA = I/Vao- 
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The average twist of the flux rope is varied between the simulations in the 
range ^ = (3-16)7r (in steps of Itt between 3 and IOtt) by varying the minor 
torus radius between a = 0.60 and a = 0.11. All these equilibria are unstable 
with respect to the helical kink mode. See [Torok, Kliem, and Titov| ([2004) for 
the calculation of the average twist, for the overall evolution of the instability 
in such equilibria, in particular for the formation of the flare current sheet and 
the helical current sheet, and for a parametric study of the growth rate. 

For equilibria of progressively higher twist, the axis of the kinked flux rope 
shows progressively higher values of the number of helical turns, with the transi- 
tion to a dominance of the double helical mode {k^ ~ 2) occurring for <l> ^ Gtt. As 
a consequence, the hehcal deformation of the rope legs leads to a progressively 
smaller separation between the legs under the rising rope apex (Figure [T|). This 
facilitates the onset of reconnection between the legs. For much higher twists 
{^ > 127r), higher axial modes (k' > 3) dominate initially, but the upward 
expansion of the flux rope leads to the dominance of the double helical pertur- 
bation in the nonlinear evolution of the instability, an effect we will address in 
a future study. 

The twist value at the transition between the regimes of single and double 
helical mode dominance depends not only upon the radial twist proflle of the rope 
as discussed in the Introduction, but also on the height proflle of the overlying 
potential fleld. This dependence will be studied in more detail in a future paper. 
Our flnding of $ > Gtt for double helical mode dominance may be compared 
with the simulation of the helical kink instability in a helmet streamer-like 
conflguration above a highly twisted flux rope in |Birn, Forbes, and Hesse|(|2006p . 
The erupting top part of their flux rope had twists in the range ~ (9-1 2) tt and 
it clearly showed the double helical shape of a /c' ~ 2 perturbation with an axis 
shape similar to the <l> = 97r snapshot in Figure [H Unfortunately, no statement 
was given whether the legs did begin to reconnect in the later phases of the 
evolution. 

2.1. Reconnection in the m = 1^ k' r^ 2 Helical Kink 

In the following we discuss one run of the series in more detail, taken to be the 
run with ^ = Qtt. This run shows a clear dominance of the double helical mode, 
which we intend to relate to the observations of the 18 April 2001 eruption. 
While this event motivated our study, we do not attempt to model the observed 
source shape, apex trajectory, and height of the crossing point very closely in this 
section. The purpose of the simulation is primarily to demonstrate that the legs 
of the rope reconnect if the twist is sufficiently high, with the apex continuing 
to rise at high speed and the region of leg-leg reconnection not participating in 
the fast rise. Quantitative comparisons between simulation and observation are 
given in Section 12.21 and in Paper I. 

Figure [2]presents snapshots of field lines in the core of the fiux rope which show 
the initial condition, the evolution of the double helical k' r^ 2 kink perturbation, 
and the approach and interaction of the rope legs. The reconnection of the outer 
fiux layers of the rope begins between the second and third snapshot pair and 
the corresponding perturbation has propagated into the core of the rope by the 



SOLA: radejection2.tex; 14 July 2010; 0:11; p. 6 



Reconnection in kinking flux rope ejection II 




Figure 2. Field lines in the core of a kink-unstable flux rope with initial average twist <l> = 97r 
(started at circles of radius r = 0.2a centered at each footpoint of the magnetic axis of the 
rope) in top and perspective views. This shows the evolution of the double helical (k' ~ 2) 
perturbation, which leads to the approach and subsequent reconnection of the flux rope legs. 
A new low- lying, less twisted flux rope is formed (bottom panels). The inner part (of size 
4 X 4 X 4) of the simulation box is displayed, including the magnetogram, Bz(x,y, 0,f). 



time of the third snapshot pair, distorting the field fines, afifiougfi tfie core stifi 
fias tfie original field line connections. Tfie final snapsfiot pair sfiows tfie system 
wfien most of tfie flux in tfie original rope fias reconnected, resulting in tfie 
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Figure 3. Left: Field lines started at the same positions as in Figure [2] at a later time, after 
the original flux rope legs have reconnected. A new flux rope is formed below the reconnection 
region, which remains at low heights (z < 1), while the top part of the original rope continues 
the upward expansion. Right: Field lines started in the top part of the rising flux at the same 
time, showing the flux rope structure and the new magnetic connections to the sources of 
the originally overlying flux. This is supplemented by some fleld lines started at the flux rope 
footpoints. The inner subvolume of size 6"^ is shown. The views differ by a rotation of 70° 
about the vertical {z) axis. 



break-up of the original rope and reformation of an only weakly kinked, still 
twisted flux rope at low heights. 

Both the rise of the rope's upper part and the reconnection at low heights 
continue subsequently. Figure [3] illustrates that the rising flux still has the 
character of a flux rope, although its average twist is reduced and although 
it becomes increasingly rooted in the ambient flux, which is obviously involved 
in the reconnection that breaks up the original rope. The downward reconnection 
outflow pushes the newly formed flux rope toward the bottom boundary, enforc- 
ing further reconnection with the ambient flux. Eventually, the newly formed 
rope breaks up as well into two flux bundles that connect each footpoint of 
the original rope to the corresponding external flux concentration of opposite 
polarity (not shown); this is analogous to the corresponding potential field. 

The rise of the apex is initially exponential, reaching O.SVao by t ^ 12rA, 
followed by an approximately linear rise (with slowly decreasing velocity), qual- 
itatively similar to the two rise profiles shown below in Figure [6] (Section 12. 2p . 
This is the typical behavior of fast CMEs and consistent with the 18 April 
2001 eruption, which showed a fast linear rise after the occulted short main 
acceleration phase. We note that the flux rope enters the height range of the 
torus instability in the course of its main acceleration. This instability, which is 
a form of the lateral kink ( Kliem and Torok, 2006D , eventually drives the flux 
rope to a full eruption. 

The isosurfaces of the current density plotted in Figure |4] show the develop- 
ment of current sheets, in addition to the kinking, interaction, and reformation 
of an approximately straight current channel in the core of the flux rope. The 
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t = 20 






t - 22 


is=m. 


^ 



Figure 4. Isosurfaces of current density, at 10% of the peak current density in each snapshot, 
showing the kinking, reconnection, and reformation of the current channel (gray) in the core 
of the flux rope which is displayed in Figures [2] and [S] In addition, the formation of the vertical 
(flare) current sheet (red) and of the lower parts of the helical current sheet (cyan) are apparent 
in panels 1 and 2. Panel 3 shows the temporary formation of a double current sheet from two 
sections of the helical current sheet between the approaching flux rope legs just before the legs 
reconnect. Upward views through the bottom of the box in panels 5 and 6 show the bottom 
edge of the squeezed current sheet (cyan) between the light-gray horizontal sections of the flux 
rope legs, the dark-cyan underside of the helical current sheet that wraps around the upper 
part of the flux rope legs (t = 15), and the current flow through the reconnected horizontal 
sections of the flux rope legs (t = 22). The inner subvolume 4^ is displayed. 



first panel displays the onset of current density steepening in the vertical (flare) 
current sheet below the rising flux rope apex (red triangular structure) and in the 
helical current sheet (two oblique cyan layers near the knees of the bent channel). 
While the vertical current sheet exists only below the rising current channel, the 
helical current sheet actually wraps around the whole channel, passing over the 
apex. Only its lower sections are visible in these isosurface plots, since the current 
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density in the helical sheet decreases rapidly with height as a consequence of 
the rapidly decreasing flux density of the ambient field. Therefore, the helical 
displacement of the upper parts of the flux rope cannot pile up much flux in 
front of the interface to the ambient field, where the helical current sheet forms. 
(See, e.g.^ Fig. 7 in Arber, Longbottom, and van der Linden, I1999[ Fig. 4 in 



IGerrard et fl/.ll200ll and Fig. 3 in [Torok, Kliem, and Titov| [2QQI1 for isosurface 



plots that show the complete helical current sheet formed around kinking fiux 
ropes which are embedded in less inhomogeneous field, and Fig. 10 in Paper I for 
field line plots of the helical current sheet. See also ILionello et aL\ 11998 for the 
formation of the helical current sheet in a representative range of kink-unstable 
fiux rope equilibria.) 

The second and third panels in Figure |4] show how the approaching legs of the 
kinking rope squeeze the two layers of the helical current sheet and the vertical 
current sheet between them into a vertical double current sheet. The fourth panel 
demonstrates that the double current sheet and the legs of the kinking current 
channel reconnect, reforming a low- lying current channel (flux rope). 

Depending on the extent of the kink and on how much the rope rises during the 
kink, one can expect the two legs of the tube to hit each other at somewhere be- 
tween right angles and an antiparallel orientation. Linton, Dahlburg, and Antiochos| 



(|2QQip explored how reconnection would happen at high twist (^ = 207r) for 
such collisions. Considering an idealized set-up with periodic boundary condi- 
tions and magnetically isolated, straight twisted flux tubes brought together by 
a prescribed slow stagnation-point flow, these authors found that antiparallel 
flux tubes of like helicity bounce off each other with hardly any reconnection 
occurring, hence they termed this case "bounce interaction." This is a natural 
consequence of the fact that the flelds at the surface of the flux tubes, which 
are dominated by the azimuthal component for such strong twist, are approxi- 
mately parallel where the tubes come into contact. They found a similar bounce 
interaction when the tubes collide at right angles with the same signs of twist 
helicity and crossing helicity (i.e., RR2: right handed tubes at a collision angle 
of 27r/8 or LL6: left handed tubes at an angle of Gtt/S, in the nomenclature of 
Linton, Dahlburg, and Antiochos[ I200ip . On the other hand, they found that 



the tubes undergo signiflcant reconnection at right angles, and in fact tunnel 
through each other, if the twist helicity and crossing helicity are of opposite 
signs (RR6 or LL2). [Lintonl (|2QQ6p then showed that this tunnel interaction can 
occur for such collisions at twists as low as $ = 97r, which equals the twist 
of the kink simulation shown in Figures [2]-|4] here. This raises the interesting 
question of whether such a tunnel reconnection could occur in the kink-induced 
reconnection simulations we present here. Unfortunately, it turns out this can 
not happen, as the kink always gives a crossing number which is unfavorable 
for tunnel reconnection. The kink instability in a left hand twisted flux tube 
induces left hand writhing, which in turn gives a left handed flux tube crossing. 
Thus the collision of flux tube legs induced by the kink instability gives the same 
sign of crossing number and twist (LL6 or RR2), and can not induce tunneling 
reconnection. 

In spite of this, we do flnd signiflcant reconnection in this kink-induced col- 
lision. This, we interpret, is due to the fact that the twist here is initially less 
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than half that of the tubes studied by Linton, Dahlburg, and Antiochos (j2QQip , 



and the kink and eruption of the tubes subsequently lower it even further. This 
allows the toroidal field to dominate over the twist field at the point where the 
tubes collide. The toroidal field is antiparallel for an antiparallel tube collision, 
and is at right angles for a right angle collision. Both orientations are favorable 
for reconnection, and so significant flux tube reconnection is generated by these 
collisions. Furthermore, the reconnection is here driven and enforced by the over- 
all helical deformation of the flux rope, which is a fast process, since the helical 
kink is an ideal MHD instability. The current channels of diameter 2a ^ 0.4 
merge completely within ^Sta (between t ^ 12.5 and t ^ 20), corresponding to 
a high reconnection inflow Alfven Mach number of ?^0.05. 

A further difference to the simulations in Linton, Dahlburg, and Antiochos 



(|2001 ) is given by the presence of the helical current sheet. This sheet and the 
current flowing in it wrap around the current channel in a helical manner, with 
the overall direction in the sheet being opposite to the direction of current flow 
in the channel. When the two layers of the helical current sheet are squeezed 
together between the approaching, nearly vertical legs of the channel, they are 
nearly antiparallel (see the third panel of Figure H]). Such a double current sheet 
is known to be susceptible to the double tearing mode. Previous two-dimensional 
simulations of this mode had sufficient spatial resolution to demonstrate that it 
creates a chain of magnetic islands in each layer, with the island growth in each 
layer facilitating the growth in the other layer. This leads to fast reconnection 
which depends only weakly on the resistivity {e.g., |Pritchett, Lee, and Drake 



fT98Ql IMatthaeus and Lamkinl [T9861 |Wang et aL\ [2QQ7|) . As in a single current 
sheet, the islands in each layer merge subsequently, driven by the coalescence in- 
stability. The merged islands — plasmoids — typically oscillate and are eventually 
ejected at high speed along the sheet if no symmetry is prescribed ( [Tajima et a/., 19 8711 

[Schumacher and Kliem, 1996[[Magara, Shibata, and Yokoyama, 1997[|Barta, Vrsfiak, and Karlicky, 2QQ8D .| 
Such acceleration of plasmoids along the sheet is amplified in the double tearing 
mode by the repelling force between neighboring plasmoids in different layers 
because the currents are oppositely directed ( [Matthaeus and Lamkin, 1986D . 
Hence, the double tearing mode has the potential to support fast reconnection, 
involving plasmoid formation and ejection, when two layers of the helical current 
sheet are squeezed between the approaching legs of a kinking ffux rope. 

Far higher spatial resolution than is available in our simulation would be 
required to resolve the details of the reconnection between the two layers of 
the helical current sheet and between the ffux rope legs in three dimensions. 
Therefore, we cannot here model the double tearing mode or other processes 
that may form plasmoids or other structures capable to conffne the accelerated 
particles in compact radiation sources, as observed in the 18 April 2001 ffare. 
Rather, we have to rely on the basic knowledge reviewed above. It must be 
kept in mind, however, that the rapid driving of the reconnection in the double 
current sheet by the approaching ffux rope legs may modify the evolution of the 
double tearing mode with respect to the quoted numerical simulations, which 
considered a double current sheet without external driving. 

The reconnection induces a peak parallel electric ffeld drop equal to the rate 



of ffux change ( [Hesse, Forbes, and Birn, 20051 , which is directly available to 
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particle acceleration, and a perpendicular electric field of order E± = O.OSVa^, 
part of which can be tapped by particles drifting in the inhomogeneous re- 
connection volume {e.g., Ambrosiano et aL\ I1988[ IKliem) I1994[ IDrake et aL\ 
I2006[ [Karlicky and Barta^ 200 7). Assuming, for example, a field strength B = 
100 Gauss and a particle density n = 10^^ cm~^, we have Vao ^ 2000 kms~^ 
and E± ~ 1 kVm~^. Assuming in addition an initial fiux rope height of 30 Mm, 
the simulation results scale to a fiux in the current channel of ~ 10^^ Mx (a 
factor < 10 lower than the fiux in the whole fiux rope), reconnected in ~ 110 s, 
which yields a parallel potential of ~ 1 GV. Both are huge values, more than 
sufficient to enable the acceleration of particles to the observed energies. 

The region of reconnection between the fiux rope legs does not participate in 
the fast rise of the fiux rope's upper part, rather it remains in the range of heights 
where it commences, z < 1, with a slight tendency to shift downwards from 
the legs' initial touching point (Figures [T]-[4]). This is in qualitative agreement 
with the observations of the 18 April 2001 eruption described in Paper I, which 
indicate that leg-leg reconnection commenced at about the time the crossing 
point of the legs had appeared at the limb, i.e., at a height slightly above 
^ O.li?0, and that the region of interaction subsequently remained near the 
limb or even retracted slightly (see also Figure [3 below) . However, the heights 
of the observed crossing point and of the commencing reconnection in the sim- 
ulation differ by a factor ~ 1.5. This value can roughly be estimated from the 
length of the polarity inversion line between the main ffux concentrations in 
AR 9415 before it rotated beyond the limb, about 10° or ~O.2i?0, which should 
be comparable to the distance, D, between the footpoints of the erupting ffux 
loop. The observationally indicated height range of leg- leg reconnection thus is 
hr ^ D/2, while the simulation above, with D = 3.3, yields h^ ^ D/3. This 
moderate disagreement can be resolved by the inclusion of a second process that 
causes the ffux to rise higher before it starts to reconnect. Candidate processes 
are an accelerated slow-rise phase, for example due to slow photospheric motions 
which lead to an amplifying inffation of the current-carrying coronal ffeld {e.g., 
ITorok and Kliemi 120031 lAulanieret a/.[|2010p . as well as any other CME driver, 
for example another instability. We consider the latter possibility in the following 
subsection. 

High twist and leg- leg reconnection do not always lead to ejective behavior. 
When the simulation is repeated with a stronger line current (as in iTorok and KliemU 
|2005) but otherwise identical parameters, a failed eruption is found. Reconnec- 
tion commences between the approaching legs of the rope, as decribed above, 
but soon thereafter it also commences in the section of the helical current sheet 
above the ffux rope apex, where the current density is now much higher due 
to stronger ffux pileup. The latter reconnection, described in ITorok and Kliemi 
(|2005j), cuts the rope at its top, thus preventing an ejective evolution. 

2.2. Helical Kink Supported by the Torus Instability 

Here we describe a ffux rope that is susceptible to the helical kink mode and to 

the torus instability from the beginning. The torus instability ( [Kliem and Torok, 2006D | 

is the process behind the well-known catastrophe of a ffux rope-arcade conffgura- 



tion ([van Tend and Kuperus, 1978 Forbes and Isenberg, 1991 Lin et a/., 1998). 
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Figure 5. Left: Field lines of a kink and torus unstable flux rope of initial twist <l> = Stt 
just before leg-leg reconnection commences, showing that the height range of leg approach is 
somewhat enlarged in comparison with the case shown in Figures [2}2] due to the additional 
upward acceleration by the torus instability. Right: The same conflguration at a later time 
showing that part of the flux in the escaping upper part of the original flux rope produces closed 
loops as a result of the reconnection between the legs, while the other part has reconnected 
with the ambient fleld. The fleld line start points lie in a circle of radius ~ (2/3)a(t) centered 
at the flux rope's apex. The inner subvolume of the box of size 6^ is shown. The views differ 
by a rotation of 70° about the z axis. 



It provides an additional upward acceleration of the flux rope, raising the height 
range of leg- leg reconnection due to the helical kink mode. In order to make 
the initial equilibrium considered in Section 12.11 also torus unstable, one has to 
modify the external poloidal field such that it falls off more rapidly with height, or 
to reduce the stabilizing influence of line-tying. Let us first describe a simulation 
that realizes the former option by reducing the distance of the external poloidal 
field sources to L = 0.25. The twist is reduced slightly to bring the time scales 
of the two instabilities closer together (so that the torus instability can influence 
the flux rope sufficiently early in the process of helical deformation), but it is 
still kept at a rather high value of ^ = Stt. The resulting evolution of the flux 
rope is very similar to the case described in Section 12. 1[ except that the upper 
part of the rising loop is now stretched upward somewhat stronger, as shown 
in the left panel of Figure [5l The reconnection commences at somewhat greater 
heights in this run, in the range z ^ 1, and the current sheets between the flux 
rope legs are squeezed in the height range up to z ^ 1.5. This remains the range 
of reconnection between the legs as the eruption evolves, as indicated by the 
shape of the field lines in the right panel of Figure [H This figure also proves the 
leg-leg reconnection by showing that part of the flux in the upper part of the 
broken rope closes with itself. The other part reconnects with the ambient field, 
forming new, less sharply defined legs. 

As in the run in Section 12. If the rise profile of the flux rope apex shows the 
two phases typical of fast CMEs: a rapid exponential rise to the peak velocity 
of nearly 0.4Va is followed by a linear further rise with very slowly decreasing 
velocity (Figure [6]). The main part of the original flux in the legs of the rope 
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Figure 6. Rise profiles of the fiuid element at the fiux rope apex for the two simulations of 
Section [Z2] that are compared to the 17 GHz images of the fiare on 18 April 2001 in Figure [3 
below (solid lines) and in Figure 9 of Paper I (dashed lines). The initial velocity perturbation 
applied at the apex is shown dotted. 



reconnects in the interval t ~(18-25)rA, i.e., shortly after the peak acceleration 
of the ejection, similar to the run in Section 12.11 and exactly as inferred for the 
18 April 2001 flare (see Paper I). Reconnection involving the newly formed legs 
and the ambient field continues subsequently at a lower rate until the run is 
terminated. 

This run reproduces the key features of the evolving shape of the microwave 
source in the 18 April 2001 fiare. Figure [7] displays, from top to bottom, the top 
section of a fiux loop, the transition to concave- upward shape, the initially slim 
ellipse, and finally the considerably expanded and presumably further rotated 
ellipse, whose bottom stays close to the limb. 

An even better agreement with the observed shape was obtained in another 
run which started from a less twisted rope (^ = Gtt), thus reducing the growth 
rate of the helical kink further and permitting a stronger infiuence of the torus 
instability. The external poloidal field in this run was chosen as in Sect ion [2Tn and 
the torus instability can grow because the line-tying is reduced by placing the 
rope at a higher position (i?=l.l,(i = 0.1), giving it a nearly semicircular initial 
shape. The comparison with the microwave images is shown in Figure 9 of Paper I 
in the same format as in Figure [7] here. The rise profile is included in Figure [6l It 
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Figure 7. Comparison of the evolving 17 GHz source in the 18 April 2001 eruption with 
the shape of the kink-unstable flux rope of initial twist ^ = Stt shown also in Figure \E\ The 
simulated flux rope is displayed at the same viewing angle as the partly occulted microwave 
source (tilted away from the observer by 26°, so that the magnetogram, Bz{x,y,0,t), is seen 
from below the simulation box). The fleld lines outline the core of the flux rope at t = 0, 12, 
17, and 24rA. 



shows the same two-phase velocity profile and the same timing of reconnection 
between the flux rope legs {t > 23rA) relative to the main acceleration phase as 
the previous two runs. 

This third run in particular brings observation and simulation in reasonable 
quantitative agreement, resolving the moderate discrepancy in the height of the 
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crossing/reconnection-onset point discussed in Section [2711 As mentioned there, 
further possibihties exist to raise the height of leg-leg reconnection onset in the 
numerical modeling; however, we refrain from further pursuing the parametric 
study of this minor aspect. 



3. Discussion 

The analysis of the microwave data of the eruptive flare on 18 April 2001 and the 
comparison with the above MHD simulations suggest that the eruption involved 
a kink-unstable flux rope that had a high twist of ^ ^ Stt. This yields a coherent 
framework to understand the inverse gamma shape of the rising main microwave 
source, as well as the triggering and ejection of superimposed compact sources 
that likely emerged from the crossing point and moved upward to the top of the 
main source. 

The inferred twist lies several times above the often quoted threshold of 2.497r 



for the helical kink instability of a line-tied flux rope (Hood and Priest, 1981 



Einaudi and Van Hoven, 1983D , so that the question arises how such a high 



value can occur. It can likely be acquired by the flux rope in the eruption 
process itself, by adding current-carrying flux to the rising and growing rope 
through reconnection under it. This addition of flux is well established theoret- 
ically {e.g., ILin and Forbes) I2000| Lin, Raymond, and van Ballegooijen] I2004p 



and observationally {e.g., |Qiu et aL\ I2004p . Resulting twist values far exceed- 
ing the above threshold have been inferred from some interplanetary CMEs 
( Gulisano et a/., 2005[ |Qiu et a/., 2007 ), but the corresponding coronal values 



remain elusive. 

It is also possible that the initial equilibrium had acquired the required twist. 
The threshold found by IHood and PriestI (I198ip and Einaud i and Van HovenI 
(|1983p is the lowest possible threshold, valid only for a uniformly twisted (Gold- 
Hoy le) flux rope of infinite radial extent. Several factors influence the threshold 
of the instability, among them the embedding of the flux rope in strong external 
current-free field, the inhomogeneity of the radial twist profile and the radius 
of the rope. In particular, it has been shown by [Batyj (|2001j) that the threshold 
can lie at twists of IOtt and higher if the radius of the rope becomes smaller 
than the pitch length of the twisting field lines. With regard to the acquired 
twist, the event on 18 April 2001 appears to be atypical, belonging to a mi- 
nority of cases with very high twist. Other candidate cases for such and even 
somewhat higher twist are described in [Vrsnak, Ruzdjak7"a nd Ro mpolt| (|199ip , 
Romano, Contarino, and Zuccarello (|2003p , and in ICho et al.\ (|2009) . 



ICho et aL\ (|2009p suggested that a kinking flux rope with reconnecting legs 
caused an eruptive X-class limb flare on 18 August 2004. They found evidence 
for the break-up of a kinked loop at the projected crossing point of the legs. 
Hot EUV-emitting plasma was produced at low heights up to the breaking 
point, presumably by reconnection, while the top part of the loop was further 
accelerated and evolved into a CME. An outward moving compact coronal hard 
X-ray source appeared at the position and time of the break-up inferred from 
EUV and Ha images; this source was seen at higher energies than the well known 
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coronal sources formed at or above the top of soft X-ray loops in some events 
(e.^., ISui and Holman( l2QQ3p . These findings are analogous to the dynamics of 
the microwave and hard X-ray sources and to the spectral extent of the hard 
X-ray source in the 18 April 2001 fiare discussed in Paper I, and they are in line 
with the simulations presented above. 

The quick reformation of a flux rope below the reconnecting legs of the original 
rope can explain the quick reformation of filaments after some eruptions, for ex- 
ample, in the 19 July 2000 event investigated by [Romano, Contarino, and Zuccarelloj 
(1201)3 ). 

Further possible cases of reconnection at or near the crossing point of rising 
kinked filaments were presented in IJi et ar](|2003p and Alexander, Liu, and Gilbertj 



(f2006) for an event on 27 May 2002 and in lLiu and Alexand er" ^2009^ for events 
on 12 June 2003 and 10 November 2004. These authors describe hard X-ray 
sources near the projected crossing point of the legs in addition to footpoint 
sources, and the latter two papers suggest leg-leg reconnection as one possible 
mechanism of the underlying particle acceleration. No break-up of the kinked 
filaments at the crossing point occurred, however, and the shapes are suggestive 
of a dominantly single helical {k^ ~ 1) kink. The kinking of twisted fiux was 
modeled for the first and third of these events in Torok and Kliemi (|2005|) and 
I Williams et aL\ (|2005|) . respectively, assuming a twist of bn in both cases, which 
led to 8i k^ ^ 1 kink. It is likely that leg-leg interaction occurred in the event 
on 27 May 2002, which showed the clearest association of the X-ray source with 
the crossing point, as a consequence of the failed nature of the eruption and the 
cut up of the filament at its top. We are preparing a corresponding simulation 
study to detail this suggestion. For the other two events, which were ejective, 
we believe that reconnection in the helical current sheet has likely played a 
role (partly similar to our suggestion here), but that the conjecture of leg- leg 
reconnection requires further study and justification. 



4. Summary and Conclusion 

MHD simulations of the helical kink instability of arched, line-tied fiux ropes 
demonstrate that the legs of the fiux rope reconnect with each other if pertur- 
bations with two helical turns, i.e., normalized axial wavenumbers k^ ~ 2, are 
dominant. This requires twists of ^ Gtt or higher. 

Leg-leg reconnection can cause a complete break-up of the original rope. 
The upper part of the original rope can evolve into a CME or stay confined 
in the corona. Its fate depends on the overlying field, as for less twisted erupting 
fiux ropes (see IRoussev et al.\ I2003| ITorok and Kliem] I2005| IKliem and Torokj 
I2006p . The lower parts can form a new fiux rope connecting the footpoints of the 
original rope. The new rope is still considerably twisted. The break-up recon- 
nection is found to involve reconnection with the ambient fiux in all simulations, 
resulting in new connections of the rising upper part of the kinked rope to the 
footpoints of originally ambient fiux. Furthermore, the downward reconnection 
outfiow pushes the reformed fiux rope towards the bottom plane, where it re- 
connects and breaks up a second time. The resulting two major fiux bundles. 
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which connect the footpoints of the original flux rope with the corresponding 
major ambient flux concentration (sunspot) of opposite polarity, represent a 
state close to the potential field of the given fiux distribution in the bottom 
plane (illustrated, e.g., in Figs. 3 and 8 of Valor i et aL\ [201 0) The simulations 
also show that the region of reconnection between the legs of the highly kinked 
fiux rope stays low in the coronal volume and that this reconnection commences 
just after the main acceleration of the fiux rope's upper part. 

Such dynamics are in overall agreement with the fiare on 18 April 2001 studied 
in Paper I, supporting the interpretation that the microwave source was formed 
by a fiux rope which was destabilized by the helical kink mode at the onset of the 
eruption's main acceleration phase and experienced leg-leg reconnection in the 
nonlinear development of the instability. The occurrence of leg-leg reconnection 
argues against a significant role for the shear field component in the writhing of 
the erupting fiux because this effect pushes the legs of the rope apart. 

A quantitative difference relative to the observations is given by the fact that 
the crossing point of the inverse gamma shaped microwave source rose to a 
height of about half the footpoint distance before the reconnection commenced, 
while the simulations yield a height of about one third the footpoint distance 
if the eruption is solely driven by the helical kink mode. This indicates that 
the initial rise, at least up to the height where the reconnection commenced, 
involved a second process in addition to the helical kink instability. The data 
do not provide clues to the nature of that process, which could have been an 
accelerated slow-rise phase or any of the triggers of fast rise discussed in current 
CME models. Simulations that include the torus instability in addition to the 
double helical kink mode, as a representative model, resolve the discrepancy. 

The approach of the fiux rope legs squeezes sections of the helical current sheet 
associated with the kink mode into a temporary double current sheet between 
them. This configuration is susceptible to the double tearing mode, which is 
known to facilitate fast reconnection and the formation of moving plasmoids 
through subsequent island coalescence. However, the limited resolution of our 
simulations does not permit us to determine whether this mode did actually 
occur. Driven by the kink instability, the double current sheet and the fiux rope 
legs reconnect essentially completely and at a high rate, with the Alfven Mach 
number of the reconnection infiow reaching ^ 0.05. 

The path of the moving compact microwave and hard X-ray sources along the 
upper part of the inverse gamma source and the coalescence of multiple blobs into 
a final compact microwave source at the top of the inverse gamma source in the 
18 April 2001 fiare are consistent with plasmoid formation and propagation in 
the helical current sheet formed by the m = 1 kink instability. The alternative 
assumption of formation by reconnection in the vertical current sheet of the 
standard fiare model appears to be incompatible with these observations (see 
Paper I). 

These findings lead us to conclude that the eruptive fiare on 18 April 2001 was 
likely caused — not exclusively but to a large part — by the helical kink instability 
of a highly twisted fiux rope. 

Acknowledgements We thank M. Shimojo for the Nobeyama Radioheli- 
ograph maps included in Figure [71 and T. G. Forbes for discussions of the 



SOLA: radejection2.tex; 14 July 2010; 0:11; p. 18 



Reconnection in kinking flux rope ejection II 



instabilities considered here. This study was supported by the DFG, by an 
STFC Rolhng Grant, by NASA Grants NNH06AD58I and NNX08AG44G, and 
by Grant 300030701 of the Grant Agency of the Czech Academy of Sciences. 
Financial support by the European Comission through the SOLAIRE network 
(MTRM-CT-2006-035484) is gratefully acknowledged. The research leading to 
these results has received funding from the European Commission's Seventh 
Framework Programme (FP7/2007-2013) under the grant agreement n 218816 
(SOTERIA project, www.soteria-space.eu). 



References 

Alexander, D., Liu, R., Gilbert, H.R.: 2006, Astrophys. J. 653, 719. 

Ambrosiano, J., Matthaeus, W.H., Goldstein, M.L., Plante, D.: 1988, J. Geophys. Res. 93, 
14383. 

Arber, T.D., Longbottom, A.W., van der Linden, R.A.M.: 1999, Astrophys. J. 517, 990. 

Aulanier, G., Torok, T., Demoulin, P., DeLuca, E.E.: 2010, Astrophys. J. 708, 314. 

Barta, M., Biichner, J., Karlicky M.: 2010, Adv. Space Res. 45, 10. 

Barta, M., Vrsfiak, B., Karlicky, M.: 2008, Astron. Astrophys. 477, 649. 

Baty, H.: 2001, Astron. Astrophys. 367, 321. 

Birn, J., Forbes, T.G., Hesse, M.: 2006, Astrophys. J. 645, 732. 

Cho, K.-S., Lee, J., Bong, S.-C., Kim, Y.-H., Joshi, B., and Park, Y.-D.: 2009, Astrophys. 
J. 703, 1. 

Drake, J.F., Swisdak, M., Che, H., Shay, M.A.: 2006, Nature 443, 553. 

Einaudi, G., Van Hoven, G.: 1983, Solar Phys. 88, 163. 

Fan, Y.: 2005, Astrophys. J. 630, 543. 

Fan, Y., Gibson, S.E.: 2003, Astrophys. J. 589, L105. 

Forbes, T.G., Isenberg, P.A.: 1991, Astrophys. J. 373, 294. 

Gerrard, C.L., Arber, T.D., Hood, A.W., Van der Linden, R.A.M.: 2001, Astron. Astro- 
phys. 373, 1089. 

Gibson, S.E., Fan, Y., Mandrini, C., Fisher, G., Demouhn, P.: 2004, Astrophys. J. 617, 600. 

Gilbert, H.R., Alexander, D., Liu, R.: 2007, Solar Phys. 245, 287. 

Green, L.M., Kliem, B., Torok, T., van Driel-Gesztelyi, L., Attrill, G.D.R.: 2007, Solar 
Phys. 246, 365. 

Gulisano, A.M., Dasso, S., Mandrini, C.H., Demouhn, P.: 2005, J. Atmos. Sol. Terr. Phys. 67, 
1761. 

Hesse, M., Forbes, T.G., Birn, J.: 2005, Astrophys. J. 631, 1227. 

Hood, A.W., Priest, E.R.: 1981 Geophys. Astrophys. Fluid Dynamics 17, 297. 

Isenberg, P.A., Forbes, T.G.: 2007, Astrophys. J. 670, 1453. 

Ji, H., Wang, H., Schmahl, E.J., Moon, Y.-J., Jiang, Y.: 2003, Astrophys. J. 595, L135. 

Karlicky, M., Barta, M.: 2007, Astron. Astrophys. 464, 735. 

Karlicky, M., Kliem, B.: 2010, Solar Phys., in press (Paper I). 

Kliem, B.: 1994, Astrophys. J. Suppl. Ser., 90, 719. 

Kliem, B., Torok, T.: 2006, Phys. Rev. Lett. 96, 255002. 

Kliem, B., Titov, V.S., Torok, T.: 2004, Astron. Astrophys. 413, L23. 

Lin, J., Forbes, T.G.: 2000, J. Geophys. Res. 105, 2375. 

Lin, J., Forbes, T.G., Isenberg, P.A., Demouhn, P.: 1998, Astrophys. J. 504, 1006. 

Lin, J., Raymond, J.C., van Ballegooijen, A. A.: 2004, Astrophys. J. 602, 422. 

Linton, M.G.: 2006, J. Geophys. Res. Ill, A12S09, doi:10.1029/2006JA011891. 

Linton, M.G., Dahlburg, R.B., Fisher, G.H., Longcope, D.W.: 1998, Astrophys. J. 507, 404. 

Linton, M.G., Dahlburg, R.B., Antiochos, S.K.: 2001, Astrophys. J. 553, 905. 

Linton, M.G., Longcope, D.W., Fisher, G.H.: 1996, Astrophys. J. 469, 954. 

Lionello, R., Velh, M., Einaudi, G., Mikic, Z.: 1998, Astrophys. J. 494, 840. 

Liu, R., Alexander, D.: 2009, Astrophys. J. 697, 999. 

Magara, T., Shibata, K., Yokoyama, T.: 1997, Astrophys. J. 487, 437. 

Matthaeus, W.H., Lamkin, S.L.: 1986, Phys. Fluids 29, 2513. 

Metcalf, T.R., De Rosa, M.L., Schrijver, C.J., et al.\ 2008, Solar Phys. 247, 269. 

Min, S., Chae, J.: 2009, Solar Phys. 258, 203. 



SOLA: radejection2.tex; 14 July 2010; 0:11; p. 19 



Kliem et al. 



Pevtsov, A. A., Canfield, R.C., McClymont, A.M: 1997, Astrophys. J. 481, 973. 

Pritchett, P.L., Lee, Y.C., Drake, J. P.: 1980, Phys. Fluids 23, 1368. 

Qiu, J., Wang, H., Cheng, C.Z., Gary, D.E.: 2004, Astrophys. J. 604, 900. 

Qiu, J., Hu, Q., Howard, T.A., Yurchyshyn, V.B.: 2007, Astrophys. J. 659, 758. 

Romano, P., Contarino, L., Zuccarello, P.: 2003, Solar Phys. 214, 313. 

Roussev, LL, Porbes, T.G., Gombosi, T.L, Sokolov, LV., DeZeeuw, D.L., Birn, J.: 2003, 

Astrophys. J. 588, L45. 
Rust, D.M.: 2003, Adv. Space Res. 32, 1895. 
Rust, D.M., Kumar, A.: 1996, Astrophys. J. 464, L199. 
Rust, D.M., LaBonte, B.J.: 2005, Astrophys. J. 622, L69. 
Sakurai, T.: 1976, Pub. Astron. Soc. Japan 28, 177. 

Schrijver, C.J., Elmore, C., Kliem, B., Torok, T., Title, A.M.: 2008a, Astrophys. J. 674, 586. 
Schrijver, C.J., De Rosa, M.L., Metcalf, T., et al: 2008b, Astrophys. J. 675, 1637. 
Schumacher, J., Kliem, B.: 1996, Phys. Plasmas 3, 4703. 
Sui, L., Holman, G.D.: 2003, Astrophys. J. 596, L251. 
Tajima, T., Sakai, J., Nakajima, H., Kosugi, T., Brunei, P., Kundu, M.R.: 1987, Astrophys. 

J. 321, 1031. 
Titov, V.S., Demouhn, P.: 1999, Astron. Astrophys. 351, 707. 
Torok, T., Kliem, B.: 2003, Astron. Astrophys. 406, 1043. 
Torok, T., Kliem, B.: 2005, Astrophys. J. 630, L97. 
Torok, T., Kliem, B., Titov, V.S.: 2004, Astron. Astrophys. 413, L27. 
Valori, G., Kliem, B., Torok, T., Titov, V.S.: 2010, Astron. Astrophys., in press, e-print 

arXiv: 1005.0254 
van Tend, W., Kuperus, M.: 1978, Solar Phys. 59, 115. 
Vrsnak, B.: 2001, J. Geophys. Res. 106, 25249. 
Vrsnak, B., Ruzdjak, V., Rompolt, B.: 1991, Solar Phys. 136, 151. 

Vrsnak, B., Ruzdjak, V., Rompolt, B., Rosa, D., Zlobec, P.: 1993, Solar Phys. 146, 147. 
Wang, Z.X., Wang, X.G., Dong, J.Q., Lei, Y.A., Long, Y.X., Mou, Z.Z., Qu, W.X.: 2007, Phys. 

Rev. Lett. 99, 185004. 
WiUiams, D.R., Harra, L.K., Brooks, D.H., Imada, S., Hansteen, V.H.: 2009, Pub. Astron. 

Soc. Japan, 61, 493. 
Williams, D.R., Torok, T., Demoulin, P., van Driel-Gesztelyi, L., Kliem, B.: 2005, Astrophys. 

J. 628, L163. 
Zhou, G. P., Wang, J. X., Zhang, J., Chen, P. P., Ji, H. S., Dere, K.: 2006, Astrophys. J. 651, 

1238. 



SOLA: radejection2.tex; 14 July 2010; 0:11; p. 20 



